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ABSTRACT /RESUME 

The paper summarizes the most important effects in 
Einsteinian gravitomagnetic fields related to propa- 
gating light rays, moving clocks and atoms, orbiting 
objects, and precessing spins. Emphasis is put onto 
the gravitational interaction of spinning objects. The 
gravitomagnetic field lines of a rotating or spinning 
object are given in analytic form. 



1. INTRODUCTION 

In the Einstein theory of gravity the gravitational field 
is described by ten potential functions g^^i, (/i, v = 
0, 1, 2, 3, gfj,i, = g^^) which at the same time are the 
metric coefficients of curved spacetime. The potential 
function goo, sometimes called gravitational redshift 
potential, is connected with the Newtonian gravita- 
tional potential, the six functions gij — 1,2,3) 
describe the geometry of the curved three-dimensional 
spaces defined by the slices of constant time t, or 
a;" = ci, where c denotes the speed of light, and finally, 
the three functions goi describe how the geometry of 
the three-dimensional slices rotates in going from one 
slice to another. For weak gravitational fields, the 
three potential functions g^i behave very much like 
the electromagnetic 3-potential A™ on account of 
which the goi are often called gravitomagnetic field 
potentials. Similarly to the generation of the elec- 
tromagnetic potential Af^ through charge currents is 
the gravitomagnetic field generated through mass cur- 
rents (momentum densities) if also by a factor of four 
more efficient because of an underlying tensor theory 
(spin-2 field theory versus spin-1 field theory of elec- 
trodynamics; implying also a sign difference). 
In this paper important effects connected with weak 
gravitomagnetic fields will be discussed. The applied 
class of gravitomagnetic fields will originate both from 
spinning and orbiting mass currents. The objects 
moving or propagating in the gravitomagnetic fields 
will be non-spinning and spinning objects (celestial 
bodies, particles, clocks, atoms, black holes, etc.) as 
well as light rays. The gravitomagnetic field lines of a 
spinning object will be given in analytic form. In all 
publications known to the author, the graphs for the 



gravitomagnetic field lines are not presented in fully 
exact form. 



2. SPINS IN MINKOWSKI SPACE 

In the weak-field limit of Einstein's theory of gravity 
the treatment of spinning objects is closely related 
to their treatment in Minkowski space. Therefore, in 
this section, properties of spinning or (rigidly) rotat- 
ing objects in Minkowski space will be discussed. The 
most important outcome of the present section is the 
relation between the canonical position variable of a 
spinning object and the various centre-of-mass defini- 
tions. 

Written in canonical variables, the total angular mo- 
mentum of a spinning object in Minkowski space takes 
the form 



J = R X P 



(1) 



where R, P, and S denote the position vector, the 
linear momentum, and the spin vector of the object, 
respectively. R and P are canonically conjugate vari- 
bles which commute with the spin vector S the compo- 
nents of which fulfil the standard angular momentum 
commutation relations. The Poincare algebra tells us 
that the centre-of-mass constant K has to take the 
form [1] 
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(2) 



where M and cP° = H are the rest mass and the 
energy of the object, respectively. The 4-momentum 
reads = (P°, P) (P^'P^ = -NPc^) and for the 
(coordinate) velocity v the relation v = cP/P'' holds. 
From the Eq. 2 the centre-of-mass coordinate results 
in the form 



using the standard definition 

K = a;°P - RP°. (4) 

In terms of the antisymmetric spin-4-tensor S^'^ 
the Eq. 4 impHes 5°* = (so-called Corinaldesi- 
Papapetrou spin supplementary condition, e.g. see 
[2]). This spin tensor we may call Si^" with S"^ = 
e^^'^Sk, where e'^*^ denotes the total antisymmetric 
Levi-Civita tensor. Then the total angular momen- 
tum 4-tensor J^" takes the form 



withX" = (X",R). 

Let us call the spin-4-tensor say, S^^'^ , if the covariant 
condition holds, S^^Ui, = (so-called Pirani spin sup- 
plementary condition, e.g. see [2]) with Uv = Pv/Mc. 
The index rf is chosen such as to indicate that the 
definition of the centre-of-mass is made in the rest 
frame. Then we get 

J^"" = X^fP" - X^fP'' + S^f", S^f^U, = 0, (6) 



TheEq. 8 also fits with electrons. The insertion of the 

spin of the electron results in the Compton wavelength 
as diameter of the area in question which is consis- 
tent with the minimum size of positive frequency wave 
functions. It is also nice to point out that even the 
radius of the ring singularity of a Kerr black hole fits 
with Eq. 8 although in Einstein's theory of gravity 
Minkowski spaces do exist only locally. 



3. GRAVITOMAGNETISM 

The line element of 4-dimensional curved spacetime 
can be decomposed into various forms (the signature 
of the metric is +2), 

ds^ = Qnvdx^dx'^ = gooc^dt^ + 2goicdtdx^ + gijdx^dx^ 

= {-goo)[-{cdt - Aidx'f + {Gij + AiAj)dx'dx^ 

= -N'^c^dt^ + gij{dx' + N'cdt){dx^ + N^cdt) , 

(9) 



with X^f 



(10f,Rrf) and 47 = e^^■'''^UJ0,i. Be- 
longing to the same reference frame, we may put 
X° = X°f = x°. The relation between R and Rrf 
is achieved by a Lorentz transformation from the rest 
frame {X^^) to the moving frame where, for centre-of- 
mass coordinates, the hat applies {X^). One finds. 



R — Rr- 



Mc PO ' 



(7) 



where Srf denotes the spin of the object in the rest 

frame, =^'"'[Skri, [3], [2]. 

Obviously, R — Rrf, i.e. the difference vector of the 
centre-of-mass positions, on the one side defined in 
the reference frame where the object is moving and on 
the other side defined in the rest frame but (Lorentz) 
transformed to the system where the object is moving, 
is orthogonal to v; thus there is no Lorentz contrac- 
tion involved. Under the assumption of positivity of 
energy density of the object in all inertial frames it 
follows that the radius of the minimum size of the 
object, orthogonal to the spin direction, is given by 
[3], [4], 



(8) 



where by definition, 
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-goo + 9ijN'N\ (11) 
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hold. The functions g^"^ define the inverse metric, 
9^^9\y — S'^- N and are called lapse and shift 
functions, respectively. 

The decomposition of the line element into the form 
where Ai appears is adapted to observers at rest in 
the given coordinate system; for those observers equal 
time means: cdt = Aidx^ , see Fig. 1. On the other 
side, the decomposition where the lapse and shift 
functions appear relates to observers at rest in "abso- 
lute spaces", i.e. in the spaces defined hy t = const., 
[5] . Those observers move in the given coordinate sys- 
tem according to dx'^ = —N^cdt, see Fig. 2. 
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Fig. 1. Shown is infinitesimal simultaneity for an 
observer located at (Einstein synchronization). 

The functions goi, Ai, and A''* are different representa- 
tions of the gravitomagnetic field (this field comprises 
the Coriolis field as well), respectively being denoted 
by g, A, and N; their relations are: goi = — goo^i = 
gikN^- Obviously, N"" = Y'^gok is valid, whereby 7'^^ 
is the inverse metric to gik, J^''gik = ^l- 




Ndt 



ds^ = - (1 1-) c^dt^ + 2ng^dtd<l) 

+ {dg^ + g'^dcp'^ + dz'^). (12) 

One easily reads off, go^ = Q.g'^/c, or, in vectorial 
notation 

g = ^nxr. (13) 

This is the well-known Coriolis field. 
In the rotating reference frame, the velocity v = 
dH/dt, with cER = dx^{t), of a particle which is at 
rest in the (global) incrtial frame connected with the 
centre of the rotating frame, is given by 

v = -Nc=-nxR. (14) 

It should be interesting for the reader to compare this 
velocity with the corresponding one in the next sec- 
tion. 

3.2 Dragging of inertial frames 

The line element of the exterior field of a rotating 
body, in non- rotating cylindrical coordinates, to lin- 
ear order in G, where G denotes the Newtonian grav- 
itational constant, reads, if the spin vector (proper 
rotation vector) S is pointing in z-direction. 
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Fig. 2. Shown are the lapse and shift functions. The 
adapted observer moves along Ndt in the spacetime. 

In the following, typical examples of Coriolis and (gen- 
uine) gravitomagnetic fields will be presented. 

3.1 The Coriolis field 

The line element of a rigidly rotating reference frame, 
in cylindrical coordinates, takes the form. 



where r"^ = g"^ + and = hold. In this case we 
find 5o<^ = -2GSg^/c^r^ or, 

2G 

g = - ^Sxr- (16) 

The velocity of a particle, as measured in the given 

coordinate system, which freely follows the action of 
the gravitational field along the azimuthal direction 
(in radial direction the particle is kept fixed), is 



V = - Nc 
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S X R + .. 



(17) 



results from the simple condition that the right side 
of Eq. 20 is an extrcmum in 3-spacc. The well-known 
Sagnac effect in optics comes from the last term on 
the right side of the Eq. 20 applied to two light rays 
which originate from a specific event and end in an- 
other one. 



This is the famous frame-dragging effect for particle 

motion; it results from the gravitomagnetic field only. 
The particles show the dragging of the inertial frames 
in the same way as small wood pieces may show the 
streamlines of flowing water. 

For more insight into the "frame-dragging" field, we 
give the orbital angular velocities of particles in circu- 
lar motion in the equatorial plane of a spinning object. 
The velocities read, e.g. see [6], 
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(18) 



The naive expectation that the particle moving in 
the dragging direction is the faster one is not cor- 
rect, rather the particle on the retrograde orbit is the 
faster (lower sign in Eq. 18). This is understandable 
from the fact that the retrograde orbit needs stronger 
centrifugal force, and thus higher velocity, to exist. 



4. SAGNAC EFFECTS 

The differential of the phase ip oi a light ray is given 
by dtp = k^dx^, where = (— w/c, k) is the wave 
4- vector with frequency u. By the aid of the disper- 
sion relation k^k^ = one can derive the integral 
representation 



V'(t,xVo,4) = 

- j ujdt + j ^(^ {Gij + AiAj)dx''dxi -|- A • rfr) , 

(19) 



where A • rfr = Aidx^. Along the light ray, where 
ds"^ = holds, the phase is constant, = 0. Hereof, 



/ ujdt= [ -{J{Gij + AiAj)dxidxi+A-dv) (20) 

Jta Jxo ^ 



follows. For time independent metric coefficients, the 
frequency w is a constant. In this case, the light path 



4.1 The Sagnac effect for atoms 

For the proper time r of an object we have c^dr'^ 
—ds^. The related action W takes the form, 



W = -McV, 



(21) 



and the quantum phase of the object is given by W/h. 
For the growth of proper time between the coordinate 
times to and t, one finds. 



T{t;to) 



to 



N^-gij{N^ + ^){m + '^) dt, 



(22) 



where = dx'^{t)/dt. The additive structure of 
cA^* -|- V* nicely reveals the motion of the reference 
frame with respect to a particle's dragged motion in 
azimuthal direction, sec section 3.2. 
The straight comparison with the change of phase in 
the light propagation case, Eq. 20, is best achieved 
through the following approximate expression. 
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(23) 



where the phase of the particle (atom, etc.) is given 

by W/h. The Sagnac effect for atoms results from Eq. 
23 in the same way as the Sagnac effect for light rays 
results from Eq. 20. By the aid of Eqs. 21 and 23 the 
clock effects originating from the gravitomagnetism 
can be deduced in a straightforward manner. 

4.2 Analogy with the electrodynamics 

The gravitational Hamiltonian of a point mass reads. 



H = Nc ^M'^c^ + Y^PiPj - cN • 



Mc'^ + (P-MNcf + M^ + (24) 

2M^ ^ ' V / 



where P = (Pi) and where the Newtonian gravita- 
tional potential $ has been introduced, goo = — 1 — 
2$/c^ + ... (expansion in powers of 1/c^). 
The analogous expressions in the electrodynamics are, 



Notice that the coordinate R is the canonical con- 
jugate to P. It relates to the particle coordinate of 
section 2 with the same name. Obviously, the total 
angular momentum 



J = L + S, 



H = c ^M^c^ + (P - -A«™)2 + 6$*^" 

= Mc^ + 4t(P - -A'^")2 + e*^'" + ... (25) 
2M c 



Obviously, the analogy, valid in the weak-field slow- 
motion limit, takes the form: eA^™ < — > MNc^, 
e$«™ < — > M$. Remember, however, the difference in 
the field- generation relations: 4Mv — > — Nc^, M — > 

— $, in the harmonic gauge, whereas ev A'^'^c, 
e correspondingly in the Lorentz gauge. 

5. SPIN DYNAMICS IN GRAVITY 

In the framework of the Einstein theory of gravity the 
gravitational interaction Hamiltonian of two spinning 
objects, linear in G and linear in each spin (in this 
approximation, the proper rotation of an extended 
object can be treated as rigid), is obtained by the 
substitution, in the momentum and Hamiltonian con- 
straints of the Einstein field equations, of the mo- 
mentum density Pa(S(x — Xq), with a ~ 1,2, through 
(Pa -h 5S0 X Va)^(x — Xq). Hereof, in the rest frame 
(vanishing total linear momentum), the spin-orbit in- 
teraction results in the form, also e.g. see [2], 
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where 



S = 
L = 

b = 



Si + 82, 
R X P, 
Si S2 



(27) 
(28) 

(29) 



and where P = Pi = — P2, R = xi — X2, R = |R| 
and for the spin-spin interaction one obtains 



G /3(Si-R)(S2-R) 



(Si • S2) 



(30) 



(31) 



is conserved in time, 

{J,Hso + Hs,s,} = 0, (32) 



dJ 

dt 



where {, } denotes the standard Poisson brackets, and 
also the absolute values of the spin vectors So (a = 
1,2) do not change with time. 



JC2 

= {SlHso + Hs,s.} = 0. (33) 

The spin-spin interaction Hamiltonian, Eq. 30, is 
identical with the spin(l)-spin(2) Hamiltonian of two 
Kerr black holes. The latter is part of the Hamilto- 
nian which describes the full spin-spin interaction of 
two Kerr black holes to linear and quadratic orders in 
G and S, respectively, [5] 
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(34) 



(35) 



It is interesting to note that the spin(a)-spin(a) inter- 
action terms (a = 1,2) are of monopole-quadrupole 
type whereas the spin(l)-spin(2) interaction is of 
dipole-dipole type. The latter is independent from 
the masses Ma, see Eq. 30. The spin-spin interaction 
Hamiltonian, Eq. 30, as well as the part of the spin- 
orbit interaction Hamiltonian, Eq. 26, which depends 
on S only result from the gravitomagnetic field, the 
other part in Eq. 26, depending on the masses, would 
also be present without gravitomagnetic field. 
For completeness we give the orbit-orbit interaction 
Hamiltonian which results from the gravitomagnetic 
field. In the rest frame, it reads, cf. [7], 
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(36) 



Again, the expression does not depend on the masses. 
5.1 The Lense-Thirring effect 



spin in the gravitomagnetic field of a spinning central 
object. It is given by 



The (proper) Lcnsc-Thirring effect is the precession 
of the orbital plane of an object moving in the grav- 
itomagnetic field of a spinning central object. By the 
aid of the orbital angular momentum vector L and the 
Runge-Lenz-Laplace-Lagrange vector [7] which is de- 
fined in the rest frame of the binary system and points 
from the centre-of-mass position to the periastron of 
the relative orbit, 



M = PxL - . , ^ . ? — . (37) 



M1 + M2R 



(M • L = 0) the precession of the orbit takes the form 



dSi 
dt 



{Si,i?sisj = f^s. xSi, (43) 



where the precessional vector reads 



G /^ (R-S,)R 



(44) 



For the gyroscopes of the GP-B mission an orbital- 
averaged spin precession of 41 mas/yr is predicted [10] 
(usually quoted as 42 mas/yr, e.g. see Fig. 5 in [11]) 
to be measured with an accuracy of 0.3%. 



(JT \ S 
— j ^ >t =< {L, Hso} >t= f^so X L , (38) 

g 

< >t =< {M, Hso} >t= J^so X M , (39) 



where <>t denotes orbital averaging and where the 
precessional frequency vector is given by 



O 1 / , (L-Seff)L \ 

S2so = ^ < >t (^Seff - 3 j , (40) 



with 



5.3 Gravitomagnetic field lines 

In the following discussion about gravitomagnetic 
field lines the object 2 with mass M2 and spin S2 
will be assumed to be much heavier than the object 
1. Then the object 1 can be treated as test object in 
the gravitational field of object 2. The shift function 
of object 2 reads 



2C 

N = — r X S2 . 



(45) 



Hereof the gravitomagnetic field strength H, analo- 
gously to the electrodynamics, follows in the form. 



Seff = S -MiMsb. 



(41) 



H = VxNc=4^rS2-3^^-|^^) . (46) 



If respectively e and a arc eccentricity and scmima- 
jor axis of the relative orbit, the averaging procedure 
yields, e.g. see [8], [7], 



^ i?3 ^* a3(l-e2)3/2- ^^^^ 



For the LAGEOS satellite, |r2so| results in 31 mas/yr 
[9]. In this case, where the index 1 may apply to the 
non-spinning satellite and the index 2 to the Earth, 
SeflF is simply given by S2. 



Obviously, 



(47) 



holds. In the exterior regime of object 2, the gravito- 
magnetic field strength allows the representation 



H 



VA 



(48) 



5.2 The Schiff effect 



with 



The Schiff effect (also called Lense-Thirring effect for 
spin, or frame-dragging effect) is the precession of a 



(49) 



As orthogonal trajectories of the A = const, lines, the 
gravitomagnctic field lines are easily obtained. The 
equation for them reads, 



(r-S2)^ 



(50) 



where Tmax IS thc maximum value of r for a given field 
line (notice: VA(r) • Vrmax(r) = 0). The gravitomag- 
nctic field strength, in the exterior regime of object 2, 
takes the form, 
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2GS2 dr{r{e),e) 
c^r-^rmax c^cos^ 



where is the angle between r and S2, i.e. r 
rS2Cos6 and where in the function r{6) the r 
to be kept constant. 



(51) 
S2 = 

max has 




Fig. 3. Shown are the gravitomagnctic field lines of 

object 2 and the spin precession of object 1 about a 
field line. The object 2 is assumed to be at rest in the 
given coordinate system, i.e. M2 Mi. 

It should be remarked that the gravitomagnctic field 
lines in Fig. 3 differ in shape from those in other 
references, [5], [6], [9], [10], [11], where [5] is closest. 

5.4 The de Sitter effect 

The de Sitter effect (also called Fokker effect, or 
geodetic precession) results from 



^SO ^ ^1 ■ 



(52) 



where 



so 



2G_ 



4Mi 



(53) 



The first term on the right side of the Eq. 53 stems 

from the gravitomagnctic field. For thc Earth-Moon 
gyroscope, i.e. the Earth-Moon binary system is re- 
garded as spinning object, it holds, 19 mas/yr, e.g. 
sec [9], and for thc GP-B mission, 6,600 mas/yr are 
expected, e.g. see [10]. In both cases Mi <C M2 ap- 
plies (mass of Earth-Moon system vs. mass of Sun, 
rcsp. mass of satellite vs. mass of Earth), so that 
only the second term on the right side of the Eq. 53 
contributes. 
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